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ABSTRACT 


A  phase-locked  loop  is  briefly  described,  together  with  the  derivation  of  the 
basic  differential  equation  which  governs  the  dynamic  behavior  of  the  loop 
during  the  pull-in  process.  The  special  case  of  the  pull-in  process  of  the 
first-order  loop  when  a  sine  wave  of  constant  frequency  is  applied  to  the 
input  of  the  loop  is  also  described.  The  relationship  between  the  frequency 
mistiming  of  the  loop,  the  initial  starting  phase  angle  of  the  input  sine  wave, 
and  the  time  required  for  the  loop  to  pull  in  is  discussed.  The  statistical 
parameters  associated  with  the  pull-in  time  is  reviewed.  In  particular, 
expressions  are  given  for  the  probability  density  function  of  the  pull-in 
time  and  the  cumulative  distribution  of  the  pull-in  time. 
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SECTION  I 


INTRODUCTION 


The  ability  of  phase-locked  loops  to  separate  narrow-band  signals  from 
wide-band  noise  has  led  to  their  increasing  use  as  FM  discriminators,  tracking 
filters,  and  synchronizing  devices,  particularly  in  space  commimications. 

It  is  possible  to  design  a  phase-locked  loop  so  that  the  loop  need  have  an 
effective  bandwidth  only  large  enough  to  pass  the  difference  between  the  input 
signal  frequency  and  the  estimate  of  this  frequency  as  generated  by  an  oscillator 
within  the  loop.  Since  this  difference  frequency  has  considerably  less  variation 
than  the  actual  input  signal,  the  loop  does  not  need  nearly  as  large  a  bandwidth 
as  would  be  needed  if  the  loop  were  merely  a  tuned  circuit  placed  between  the 
system  input  and  the  output,  which  would  have  to  pass  all  frequencies  over 
which  the  input  signal  varied. 

Since  the  bandwidth  of  the  phase-locked  loop  can  be  made  much  smaller 
than  that  of  a  comparable  nontracking  filter,  the  amount  of  noise  reaching  the 
output  is  considerably  smaller,  and  the  loop  develops  a  greater  resistance  to 
noise  at  the  input. 

However,  improvement  in  one  aspect  of  a  system  usually  results  in  a 
sacrifice  in  another,  and  phase-locked  loops  are  no  exception.  The  reduction 
in  the  effective  noise  bandwidth  results  in  an  increase  in  the  pull-in  time  of  the 
loop,  that  is  an  increase  in  the  time  between  the  appearance  of  a  signal  at  the 
input  and  satisfactory  operation  of  the  loop. 

The  pull-in  time  can  usually  be  neglected  for  systems  having  long  periods 
of  continuous  communication  time  available,  as  would  the  case  for  a  stationary 
ground-based  system  with  direct  line  of  sight.  However,  instances  do  arise 
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where  the  pull-in  time  of  the  loop  may  be  a  significant  fraction  of  the  commimica- 
tion  time  available.  Two  examples  that  come  to  mind  are  a  meteor- burst  scatter 
system,  where  the  ionized  trails  only  last  for  a  few  seconds,  and  a  military 
satellite  commimication  system  subject  to  pulse  jamming,  where  a  straight¬ 
forward  defense  against  the  jammer  is  achieved  by  gating  the  satellite  receiver 
on  and  off  in  a  pseudo-random  fashion. 

With  these  comments  in  mind,  it  is  apparent  that  the  pull-in  time  of  a 
phase-locked  loop  may  be  a  significant  factor  in  the  design  of  certain  communi¬ 
cation  systems.  To  illustrate  the  problem  of  pull-in  time,  the  example  of  the 
simplest  type  of  phase-locked  loop  is  chosen,  and  the  time  required  for  the 
loop  to  pull-in  tvhen  an  input  signal  of  constant  frequency  is  applied  to  the  input 
is  derived. 
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SECTION  n 


DESCRIPTION  OF  PHASE-LOCKED  LOOP 


The  basic  elements  of  a  typical  phase-locked  loop  are  shown  in  Fig.  1 
and  consist  of:  a  voltage  controlled  oscillator  (VCO)  with  a  nominal  frequency 
equal  to  that  of  the  input  to  the  loop;  a  phase  detector  which  compares  the  phase 
of  the  output  of  the  oscillator  with  that  of  the  input  to  the  loop;  an  amplifier 
which  amplifies  the  output  voltage  of  the  phase  detector;  and  a  low-pass  filter 
which  filters  the  output  voltage  of  the  amplifier  before  it  is  applied  to  the  VCO. 

Briefly,  the  phase-locked  loop  operates  as  follows:  The  phase  detector 
beats  the  signal  input  and  the  VCO  output  together,  giving  a  low-frequency 
output  proportional  to  the  sine  of  the  phase  difference  between  the  two  signals, 
together  with  a  high-frequency  component  located  at  the  sum  frequency  of  the 
two  inputs.  The  low-pass  filter  accepts  only  the  low-frequency  term,  which 
is  applied  as  a  control  voltage  to  the  VCO,  forcing  the  output  phase  of  the  VCO 
to  follow  the  input  signal  phase.  Figures  2  and  3  show  the  control  characteristics 
of  the  phase  detector  and  the  voltage-controlled  oscillator.  The  phase  detector 
generates  a  voltage  (t)  which  is  proportional  to  both  the  amplitudes  of  the 
two  inputs  signals  and  to  the  sine  of  their  phase  difference,  the  constant  of 
proportionality  being  k^.  (K^  has  dimensions  of  (volts  x  radians)  ^. )  The 

U  U 

voltage -controlled  oscillator  shifts  its  output  frequency  by  an  amount  propor¬ 
tional  to  the  voltage  applied  to  its  control  input,  the  constant  of  proportionality 
being  (radians/sec /volt).  The  point  about  which  the  frequency  of  the 
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VCO  shifts  is  cug  ,  usually  referred  to  as  the  natural  or  free-running  frequency 
of  the  VCO.  Figure  3  indicates  that  co^  is  the  output  frequency  of  the  VCO 
when  the  control  voltage  is  zero. 

The  differential  equation  governing  the  d5Tiamic  behavior  of  the  loop  is 
derived  below.  The  input  signal  to  the  loop  is  assumed  to  be  VT  E  sin  (wq  +  (t)), 

where  E  is  the  rms  amplitude  and  (t)  represents  the  variation  of  its  phase 
as  a  fimction  of  time.  The  output  of  the  voltage-controlled  oscillator  is 
assumed  to  be  VT  C  cos(a)^  +  62  (t)),  where  C  is  the  rms  amplitude  and 
represents  the  variation  of  its  phase  as  a  fimction  of  time. 

A  natural  question  at  this  point  is  why  are  the  loop-input  frequency  and 
the  VCO  frequency  expressed  in  terms  of  a>^  plus  a  time-varying  phase  angle 
rather  than  in  terms  of  two  different  frequencies,  for  instance,  <01  and  (02. 
When  the  loop  input  is  frequency  modulated,  both  the  input  and  VCO  frequencies 
change  as  a  function  of  time  during  and  after  the  condition  of  lock-on.  If  the 
input  frequency  is  held  constant,  the  VCO  frequency  changes  as  a  function  of 
time  during  the  lock-on  process,  and  then  remains  constant.  Therefore,  both 
frequencies  generally  vary  as  a  function  of  time.  Hence,  it  is  very  convenient 
(in  a  mathematical  sense)  to  express  both  frequencies  in  terms  of  a  frequency 
which  is  fundamental  to  the  system  and  which  always  remains  constant.  The 
natural  frequency  of  the  VCO,  is  just  such  a  frequency.  For  example,  the 
input  to  the  loop  was  written  as 

V'2  E  sin  (t)]  (1) 

Suppose  the  loop  input  is  a  sinusoid  of  constant  frequency  given  by 

E  sin  +  a) 
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We  could  just  as  correctly  write  it  as 


V  2  E  sin  +  (<i>j  —  (u^)  t  +  a] 


(2) 


Which  when  compared  with  Eq,  (1)  gives 

6^  (t)  =  t  +  a 

Therefore,  there  is  nothing  lost  by  expressing  signals  of  different  frequencies 
all  in  terms  of  the  same  frequency,  provided  the  necessary  adjustment  is  made 
in  the  time-varying  phase  angle. 

Continuing  the  development  of  the  differential  eqxiation,  the  output  of  the 
phase  detector  is  given  by 

Cj  (t)  =  \/~2  E  sin  +  0j(t)  ]  \/~T  C  cos  +  02(')1 

(3) 

=  EC  1  sin  +  6^(t)  +  02(^)i  +  sin  [0j(t)  - 

The  term  2  in  Eq,  (3)  can  be  neglected,  since  neither  the  low-pass 
filter  nor  the  VCO  will  respond  to  it,  provided  is  reasonable  large. 

The  output  of  the  filter  is  then  given  by 

e^.(t)  =  Ki  F(S)  EC  sin  [0j(t)  -  02(t)]  ^4) 

where  F(S)  represents,  in  operational  notation,  the  effect  of  the  filter  on  the 
signal  sin  (t)  -  02  (t)]  • 
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Let 


0l(t)-02W  =  ^(0 


(5) 


where  0(t)  is  the  instantaneous  phase  difference  between  the  input  signal  and 
the  output  of  the  VCO. 

Equation  (4)  then  becomes 


e^(t)  =  Kj  F(S)  K,9  EC  sin0(t) 

The  output  frequency  of  the  VCO  is  given  by 

oJ^  +  02(^)  =  +  ^vco 

Substituting  for  e^.  (t)  from  Eq.  (6)  gives 

^2^')  =  ^vco  ^0  EC  F(S)  sin0(t) 


Differentiating  Eq.  (5)  with  respect  to  time  gives 


^2(0  =  0i(t)-d(t) 

Substituting  in  Eq.  (7)  for  02  (t)  gives 

0l(t)  -  0(t)  =  Kj  EC  F(S) 


(6) 


(7) 


(8) 


(9) 
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or 


0(t)  =  0i(t)  -  Kj  Kg  EC  F(S)  sin0(t) 


(10) 


Let 


Kfl  K,  EC 


Equation  (10)  may  now  be  rewritten  as 


(9(t)  =  (9i(t)  -  KF(S)  sin(9(t) 


(11) 


Equation  (11)  is  the  basic  differential  governing  the  behavior  of  the  phase-locked 
loop  where 


d(t)  =  instantaneous  frequency  difference  between  the  input  to  the 
loop  and  the  output  of  the  VCO. 


0^  (t)  =  instantaneous  frequency  difference  between  the  input  to  the 


loop  and  the  free-running  frequency  of  the  VCO. 

0(t)  =  instantaneous  phase  difference  between  the  input  to  the  loop 

and  the  output  of  the  VCO. 

The  terms  "instantaneous  phase  difference"  and  "instantaneous  frequency 
difference"  are  discussed  in  Appendix  V. 
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SECTION  m 


CONSIDERATION  OF  FIRST-ORDER  PHASE- LOCKED  LOOP 

The  simplest  type  of  phase-locked  loop  is  that  shown  in  Fig,  1,  except 
that  the  low-pass  filter  is  omitted,  a  direct  connection  being  made  between 
the  output  of  the  amplifier  and  the  input  to  the  VCO  [F(S)  =  1]  .  Althou^  the 
low-pass  filter  is  not  present  in  the  feed-back  loop,  it  is  assumed  that  the  VCO 
will  not  respond  to  high-frequency  components  applied  to  its  input  (that  is,  the 
term  in  2  resulting  from  product  detection). 

We  now  consider  the  case  where  a  sine  wave  of  constant  frequency  cai  is 
applied  to  the  loop. 

For  the  case  discussed  above,  Eq.  (11)  takes  the  particularly  simple 
form  given  by 


0(t)  =  Aco  —  K  sin0(t) 


(12) 


where 

0(t) 

instantaneous  phase  difference  between  the  input  to  the 

loop  and  the  output  of  the  VCO. 

m 

=  instantaneous  frequency  difference  between  the  input 

to  the  loop  and  the  output  of  the  VCO. 

Ao)  =  CJj  — 

=  frequency  difference  between  the  input  to  the  loop  and 

the  free-rimning  frequency  of  the  VCO. 

K 

=  open- loop  gain  of  the  loop. 

9 


Equation  (12)  is  a  differential  equation  of  first  order,  so  that  the  phase- 
locked  loop  shown  in  Fig.  1  (with  the  low-pass  filter  omitted)  is  usually  referred 
to  as  a  first-order  loop.  Basically,  Eq.  (12)  states  that  when  a  sine  wave  of 
constant  frequency  is  applied  to  the  first-order  loop,  the  sum  of  the  instantaneous 
frequency  difference,  plus  K  times  the  sine  of  the  instantaneous  phase  differ¬ 
ence,  is  always  constant,  that  is 

0(t)  +  K  sin0(t)  =  Ao)  (13) 

In  general,  there  will  be  both  a  frequency  difference  and  an  initial  phase 
difference  between  the  input  to  the  loop  and  the  free-running  frequency  of  the 
VCO.  The  process  of  pulling  the  VCO  into  frequency  alignment  with  the  input 
to  the  loop  will,  therefore,  require  finite  time.  Further  on  in  this  report,  the 
pull-in  time  is  shown  to  be  a  function  of  both  the  initial  frequency  detuning 
(wj  -  <Uq)  and  the  initial  phase  angle  0- . 

A  plot  of  Eq.  (12)  helps  in  understanding  the  physical  process  by  which 
the  phase-locked  loop  pulls  into  frequency  alignment  with  the  input;  it  is  shown 
plotted  in  Fig.  4. 

The  first  point  to  note  is  that  the  system  represented  by  Eq.  (12)  is  in 
equilibrium  when  the  instantaneous  frequency  0(t)  is  zero. 

This  occurs  when 

Ao)  -  K  sin0  =  0  (14) 

or 

A<u 

sin0  =  - 

K 
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*  I  |iHll||ifP|||)|  ffiff  ^wy 


1 


Fig.  4.  Plot  of  0  Versus  0  for  a  First-Order  Phase-Locked  Loop 


Let 


sin 


-1 


A&i 

K 


=  0„ 


Then  the  system  defined  in  Eq.  (12)  is  in  equilibrium  for  values  of  0 
given  by 


0  =  nrr  +  (—1)*'  0^ 


(15) 
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where  n  is  an  integer  including  zero. 

The  values  of  d  given  by  Eq.  (15)  obviously  occur  at  the  intersections 
of  the  curve  with  the  0(t)  =  o  axis. 

Which  of  the  values  of  6  given  by  Eq.  (15)  represent  positions  of  stable 
equilibrium?  For  6  positive  the  rate  of  change  of  6  with  time  is  positive, 
implying  that  6  is  increasing.  On  the  other  hand,  6  negative  implies  that  6 
is  decreasing  with  time. 

These  two  conditions  indicate  that  for  any  phase  position  for  which  the 
curve  is  above  the  6  =  0  axis  (Fig.  4),  6  moves  to  the  right  (6  increasing), 
while  for  phase  positions  below  the  6=0  axis,  6  moves  indicate  the  direc¬ 
tions  in  which  6  moves  with  time.  These  conditions  indicate  that  values  of  6 
in  Eq.  (15)  for  which  n  is  even  are  points  of  stable  equilibrium,  while  those 
for  which  n  is  odd  are  points  of  unstable  equilibrium.  Since  the  curve  drawn 
in  Fig.  4  is  periodic,  it  is  sufficient  to  restrict  our  attention  to  any  convenient 
27t  radians. 

If  the  mistuning  Aco  is  increased,  the  phase  trajectory  shown  in  Fig.  4 
moves  upwards  with  respect  to  the  6  =  0  axis  until  the  two  values  of  6  (stable 
condition  and  unstable  condition)  coincide.  For  this  value  of  6,  say  ^cric  * 

A(o  _  j  Aco 

-  =  n  —  sin  - 

K  K 

Aco  77 
K  "  2 

Aco 
K 


or 


-1 


cnc 


-1 
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Hence,  for  the  case  where  the  mistuning  A<u  equals  the  open-loop  gain  K,  the 
system  reaches  a  position  of  unstable  equilibrium  where  any  slight  disturbance 
of  the  phase  due  to  noise  within  the  loop  may  cause  the  phase  to  recycle  to  the 
next  point  of  unstable  equilibrium. 

The  above  considerations  lead  to  the  conclusion  that  the  mistuning 
A<u  =  (<uj  -  <Uq)  should  be  less  than  K  rads/sec  so  that  the  first-order  phase- 
locked  loop  settles  down  to  a  steady  condition;  i.  e. , 

|A<u|  =  |<i>j  —  tUjjl  <  K  radians/sec  (16) 

For  the  condition  given  by  Eq.  (16),  the  steady-state  phase  error  0^  is  given 

by 


.  A&j 
6  =  sin  - 

°  K 
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SECTION  IV 


PULL-IN  TIME  FOR  LOOP 


The  time  required  for  the  loop  to  pull-in  is  given  by  the  time  required  for 
the  phase  difference  6  to  traverse  the  phase  trajectory  from  some  initial  value 
to  its  steady-state  value  6^  .  The  time  required  to  reach  6^  is  always 
infinite  and  independent  of  ,  since  as  6  approaches  ,  6  changes  very 
slowly  (i.  e, ,  d  approaches  zero,  see  Fig.  4).  This  point  is  discussed  fully 
in  Appendix  1. 

It  is  more  realistic  to  say  that  pull-in  time  is  the  time  required  for  6  to 
change  from  some  initial  value  to  (Oq  -  e)  or  (0^  +  e)  depending  upon 
lying  to  the  left  or  right  of  6^  (see  Fig.  4).  Appendix  I  shows  that  the  pull-in 
time  T  is  given  by 


for  (6^  +  ()<9-<(Tr-  9^)  . 
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If  the  initial  value  of  lies  in  the  interval  (0^  -  t)  <  <  0^  +  t,  the  pull-in 

time  is,  of  course,  zero,  since  the  phase  angle  is  already  within  the  region 

defining  the  condition  of  lock-in.  Stated  formally,  T  =  0  ,  for  (d^  -  e)  < 

• 

Examination  of  Eqs.  (17)  and  (18)  indicates  that  the  pull-in  time  T  is 
infinite  when 

secd^  =  oc  or  =  —  ;  this  corresponds  to  the  case  already 

discussed,  in  which  the  detuning  Aw  is  increased  until  it  equals 
the  open-loop  gain  K . 

K  =  0  ;  this  is  a  trivial  solution  corresponding  to  the  open-loop 
gain  K  =  0. 

€ 

sin  —  =  0  or  «  =  0  ;  this  has  also  been  discussed  previously;  it 

takes  infinite  time  for  the  loop  to  pull-in  to  its  steady-state  phase 
error  0^  . 

+  ^i  +  ^o  X. 

-  =  0  ,  or  -  =  —  (  n  odd).  This  occurs  when 

2  2  2 

0.  =  (nn  -  6^)  and  corresponds  to  the  initial  phase  angle  starting 
at  the  points  of  unstable  equilibrium  given  by  Eq.  (15)  for  n  odd; 
i.e.,  =  (^  +  ^i  ^ 

zero  when 

0-  =  (0Q  -  e)  for  Eq.  (17)  and  0-  =  (d^  +  t)  for  Eq.  (18). 

The  actual  definition  of  pull-in  time  is  somewhat  arbitrary,  since  it 
depends  on  a  satisfactory  choice  being  made  for  the  value  of  « .  In 
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this  report,  the  time  taken  for  the  loop  to  pull-in  to  within  5  degrees 
of  its  steady  state  value  (e  =  5  degrees)  will  be  called  the  pull-in 
time. 

The  expressions  for  the  pull-in  time  T  given  by  Eqs.  (17)  and  (18)  are 
rather  involved;  however,  they  show  that  T  is 

inversely  proportional  to  the  open-loop  gain  K,  (For  the  first-order  loop 

K 

discussed  in  this  report,  the  noise  bandwidth  is  equal  to  —  cps  ,  so 
that  the  pull-in  time  is  also  inversely  proportional  to  the  noise  bandwidth. ) 
a  nonlinear  function  of  the  initial  phase  angle 

a  function  of  the  steady-state  phase  angle  6^  ,  which  in  turn  is  a  function 
of  the  ratio: 


frequency  mistuning  Acu  /  A(u\ 

- : -  =  -  li-CM  =  - 1 

open-loop  gain  ^  \  ^  / 

Equations  (17)  and  (18)  can  be  rearranged  so  that  the  product  T  •  K  is 
■  then  a  function  of  the  steady-state  phase  angle  6^  and  the  initial  phase  angle  0.^ . 

Figure  5  shows  T  •  K  plotted  against  the  initial  phase  angle  0.  , 
for  dg  =  0  ,  30  and  60  degrees  corresponding  to 


AdJ 

-  =  0,0.5  and  0.866. 

K 


(19) 


Figure  5  also  shows  that  certain  characteristics  regarding  the  pull-in 
time  T  now  become  apparent. 
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T 


T 


di 


Fig.  5.  Plot  of  Pull-In  Time  x  Op€n-Lx)op  Gain  (T  x  K)  Versus  Initial  Phase 
Angle  0^  for  0^  =0°,  30°  and  60° 
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T  starts  to  increase  rapidly  for  initial  phase  angles  approaching  -(n  +  6^) 
and  (n  -  0  )  for  all  three  values  of  e  ■  The  values  6:  =  -(n  +  6  )  and 

o  o  o 

=  in  -  0^)  are  actually  asymptotes  corresponding  to  infinite  pull-in 
time. 

The  two  branches  of  the  curve  for  =  0  are  symmetrical  about  the 
line  6i  =  0^  =  0  ;  as  the  mistuning  increases  (0^  increasing),  the  two 
branches  become  asymmetrical  about  the  line  •  There  is  also 

a  general  increase  in  the  pull-in  time  as  6^  increases. 
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SECTION  V 


SUMMARY  OF  PROBABILITY  CONSIDERATIONS  ASSOCIATED  WITH 
PULI^m  TIME  T  FOR  FIRST-ORDER  PHASE-LOCKED  LOOP 


The  results  obtained  in  Appendixes  I,  II,  in,  and  IV  regarding  the  statisti¬ 
cal  nature  of  the  pull-in  time  are  discussed  in  this  section. 

In  Appendix  I  (Eqs,  (46) and  (50))  it  is  shown  that  the  pull-in  time  T  is 
given  by 


(20) 


for  -  (;r  +  d^)  <  6-  <  (6^  -  e); 
and  by 


(21) 


for  6^  +  (  <  d-^  <  (n  -  d^). 
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Let 


KT  =  X 


cos 

i)] 

€ 

sin  — 

2 

and 


€ 

sin  ““ 
2 


Equations  (20)  and  (21)  become 


X  ==  secd^ 


:  B  +  In 


sin 


cos 


for  -(tt  +  6^)  <  0-  <  (0^  -  £)  ; 


and 


(23) 


for  (G^  +  f)  <  (9j  <  (;7  -  e^) , 

The  probability  density  for  x  (derived  in  Appendix  II)  is  given  by 


p(X)  =  — 


77 


(X  =  0) 


(24) 


and  by 


cos^  G, 


p(X)  = 


[= 


2  77  I  cosh(cos0^X— B)— sin0^ 


(X  >  0) 


cosh(cos0  X— B— )+sin0 


] 


(25) 


The  probability  density,  as  given  by  Eq.  (25),  is  shown  plotted  in  Fig.  6 

/^co  \ 

for  G^  =  0,  30  and  60  degrees  (—  =  0,  0.  5  and  0.  866  j  .  The  value  of  ( 
was  again  taken  as  f  =  5  degrees. 
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Fig.  6.  Probability  Density  Function  for  the  Product  (T  x  K) 


The  mean  of  X,  (X) ,  (Eqs.  (98)  and  (99),  Appendix  HI)  is  given  by 


X 


1 

-  secd^ 

2n  ° 


B(77  +  26^  —  e)  +  (77  —  26^  —  0  + 


2 


In  (cosA)dA  —  4 


In  (sinA)dA 


(26) 
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If  (  is  small,  the  last  two  terms  in  Eq.  (26)  may  be  ignored  and  the  equation 


reduced  to 


1  1 
X  «  —  secfl  [B(n+26J+B‘-(n-2d)] 


(27) 


For  f 


small  B  =  =  In  (cos0  cot  IV  SO  that 

\  V 

/2  cosfl  \ 


X  =  seed 


(28) 


for  e  small. 

Figure  7  shows  how  x  varies  as  a  function  of  the  steady-state  phase 

Aco  _ ,  — 

angle  6^,  which  is  related  to  Aco  and  TK  by  -  =  sin  ^  0^  .  X  has  been 

plotted  against  6^  for  c  =  5  degrees  and  0.  5  degrees,  and  indicates  how  the 
mean  value  of  the  pull-in  time  increases  as  the  loop  pulls  in  closer  to  the 
steady-state  phase  angle. 

The  cumulative  distribution  for  X  is  derived  in  Appendix  IV  and  is 
given  by 


P(X  <  X^) 


Figure  8  shows  the  cumulative  distribution  of  X  for  6^  =  0,  30 
60  degrees.  The  value  of  c  is  again  taken  as  5  degrees. 


P.  Bratt 


and 
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24 


Fig.  8.  Cumulative  Distribution  for  X,  ^  =  5° 
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APPENDIX  I 


DERIVATION  OF  THE  FORMULA  RELATING 
PULI^IN  TIME  T  TO  THE  INITIAL  PHASE  ANGLE  6-^ 

The  time  required  for  the  phase  difference  in  the  loop  0  to  traverse 
the  phase  trajectory  shown  in  Fig.  4  between  an  initial  value  0-^  and  a  final 
value  9{  is  obtained  by  integrating  Eq.  (12)  and  evaluating  for  0]  and  9( . 
Equation  (12)  may  be  rewritten  as: 


K  dt 


sin^ 


K 


or 


K  dt 


where 


a 


(30) 


K 


Integrating  Eq.  (30)  gives 


K  dt 
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or 


yj  \  -  a? 


In 


1  _  (l  + 


a  tan 


-  (l  -  Vl  -a^) 


=  KT 


0: 


(31) 


Let 


1 

—  In 
A 


a  tan  —  —  (1  +  A) 
2 


a  tan  —  —  (1  —  A  ) 
2 


KT 


Substituting  the  limits  in  Eq.  (32)  gives 


(32) 


1 

A 


a  tan  —  —  (1  +  A) 
2 


a  tan  —  —  (1  —  A) 
2 


—  In 


a  tan  —  —  (1  4-  A) 
2 


a  tan  —  —  (1  —  A) 
2 


=  KT 


(33) 


Before  going  further  with  this  equation,  consideration  should  be  given  to  the 
range  of  values  of  0^  and  for  which  evaluation  of  Eq.  (33)  is  of  interest. 
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Basically,  we  reqxiire  the  time  taken  for  0  to  change  from  an  arbitrary 
initial  value  0-^  to  its  final  value  Of  =  sin“^  a  =  0^  ,  which  is  the  steady- 
state  phase  error,  as  explained  previously.  Therefore  Of  in  Eq.  (33)  is  simply 
0^ ,  the  steady-state  phase  error. 

We  now  consider  the  range  of  initial  values  of  Of  to  be  used  in  evaluating 
Eq.  (33).  Figure  4  shows  that  0.  falls  into  two  convenient  ranges.  They  are 

-(tt  +  0„)  <0f<0^  (34) 


and 


0o  <  0i  <  ;r  -  0o  (35) 

Returning  now  to  the  evaluation  of  Eq.  (33),  if  the  value  0^  is  substituted 
for  Of  in  the  first  term  of  the  left-hand  side  of  the  equation  the  denominator 
becomes  zero  and  the  term  itself  becomes  infinite.  This  arises  because 

1  -  A 

tan  —  =  -  (sinO^  -  a  and  cosO  =  A) 

o  o 

a 

d0 

This  is  to  be  expected,  since  as  0  approaches  0^,  —  tends  to  zero; 

dt 

e.  g. ,  the  phase  is  changing  very  slowly,  requiring  infinite  time  to  reach  its 
steady-state  value  of  0^.  Therefore,  the  time  required  for  the  phase-locked 
loop  to  pull  into  its  steady-state  value  becomes  infinite,  independent  of  the 
initial  value  of  the  phase.  To  avoid  this  difficulty,  the  time  required  for  the 
phase  0  to  pull-in  to  within  e  of  the  steady-state  value  0^  from  an  arbitrary 
value  of  Of  is  determined  (see  Fig.  4). 
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Equation  (33)  is  evaluated  below  in  two  separate  sections  according  to  the 
two  possible  ranges  of  as  given  previously,  i.  e. , 

-{tt  +  ^o)<  <  (^o  ” 


(0O  +e)<d-<(7T-  d^) 


PULL-IN  TIME  FOR  -{n  +  B^)  <  d-^<  -  e 

The  first  term  in  Eq.  (33)  is  L  where 


L 


In 


a  tan  —  —  (1  +  A) 
2 


a  tan  —  —  (1  —  A) 
2 


(36) 


where 


Expanding 


(0o  -  0 

tan  —  =  tan  - 

2  2 


e 

—  tan  — 
2 


1  +  tan 


2 


1 

.  tan  — 

2 


30 


Now 


(1  -  A) 

tan  —  =  - 

2  a 


and  let 


tan  —  =  a 
2 


(Oo  -  0 

tan  —  =  tan  - 

.  .  2  2 


(1  -  A)  -  a  a 
a  4-  (1  —  A)a 


(37) 


Substituting  for  tan  — 


in  Eq.  (36)  gives,  after  some  simplification, 


L  = 


In 


( 


«  (A  +  a  a) 
a  (1  -  A) 


) 


(38) 


We  now  consider  the  second  term  in  Eq.  (33),  M,  where 


M 


In 


a  tan  —  —  (1  A) 
2 

a  tan  —  —  (1  —  A) 
2 


(39) 
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When  considering  the  initial  values  of  0 ,  ,  it  is  convenient  to  use  the  point 

of  unstable  equilibrium  at  -  {n  +  6^)  as  a  reference  point  in  the  same  way  that 
the  final  values  of  0  were  measured  with  reference  to  the  stable  point  0  =  0^ 

Referring  to  Fig.  4  and  measuring  from  the  point  of  imstable  equilibrium, 
it  is  seen  that 


0j.  =  (77  +  0^)  +  8 


=  —  n  —  (0  —  8) 


tan  — 
2 


tan 


(00  -  5) 


cot 


0  ^ 

o  O 

1  +  tan  —  tan  — 

2  2 

r 

tan  —  —  tan  — 

2  2 


(1  -  A)  S 

Recalling  that  tan  —  =  -  ^  and  letting  tan  —  =  ^3  , 

2  a  2 


0: 


1  + 


1  -  A 


1  -  A 


-/3 


a  +  (1  -  A)  /3 
(1  -A)-a/3 


(40) 


Substituting  in  Eq.  (39)  for  tan  —  ,  we  obtain,  after  simplification. 


M  =  In 


a/3 


t  (1  -  A)  (A  +  a  /3) 


(41) 
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Using  the  values  obtained  for  L  and  M  from  Eq.  (38)  and  (41)  and  substituting 
in  Eq.  (33)  gives 


Recalling  that 


Recalling  that 


€ 

a  *=  tan  — 
2 


and 


j8  =  tan 


8 

2 


Eq.  (42)  becomes 


KT  = 


1 

A 

1 

KA 


(42) 


(43) 


A  =  cos  6^  and  a 


sin  6^,  Eq.  (43)  may  be  written  as 


(44) 
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Also 


0[  =  -(ff  +  0^)  +  S 

for  -  (!7  +  0^)  <  -  f,  so  that 

S  =  d.^  +  n  +  d^  (45) 

Equation  (44)  may  now  be  rewritten  by  substituting  the  dependent  variable  0-^ 
for  S  . 


Equation  (44)  now  becomes 


1 

T  =  _ 

K  cos  0 

o 


In 


€ 

sin  _ 

2 


(46) 


PULL-IN  TIME  FOR  (d^  +  ()  <  0-^  <  (tt  -  6^) 

Using  the  same  approach  used  in  evaluating  T  for  the  range 
-  (n  +  d^)  <  <  (0Q  -  f),  but  using  the  point  d  =  {n  -  6^)  as  a  reference 

for  the  initial  value  of  0  (see  Fig.  4),  it  can  be  shown  that 


(47) 
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which  reduces  to 


T  = 


(48) 


We  can  rewrite  Eq.  (48)  with  6.^  as  the  dependent  variable  since 

e.  =  n-e^-s , 

for  (0Q  +  f)  <  <  (77  -  d^),  or 

5  =  {.n-d^-6-^) 


Therefore  Bq.  (48)  becomes 


1 

T  =  - 

K  cos  d 


(49) 


(50) 
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APPENDIX  II 


PROBABILITY  DENSITY  FOR  THE  PULL-IN  TIME 
FOR  A  FIRST-ORDER  PHASE -LOCKED  LOOP 


In  Appendix  I  it  was  shown  that  the  pull-in  time  for  the  first-order 
phase-locked  loop  is  given  by  (see  Eqs.  (43)  and  (47)). 


KT  =  -  In 

A 


Cf  ^  \  ^  \1 

■  A  cot  —  -I-  a  I  |A  cot  —  +  a  II 

|V  2  y  \  2/1 


for  +  Qq)  <  ■> 

and 


1 

KT  =  — 
A 


(51) 


(52) 


for  (0Q  +  €)  <  . 

Equations  (46)  and  (50)  gave  the  pull-in  time  with  the  initial  phase  angle 
0.  as  the  dependent  variable.  In  this  and  later  calculations  on  probabilities, 
it  was  found  that  the  mathematics  was  greatly  simplified  by  using  5  as  the 
dependent  variable.  It  is  then  quite  simple  to  substitute  for  5  in  terms  of  9-^ 
at  the  end  of  the  calculations.  For  convenience  we  let  KT  =  x  in  Eqs.  (51) 
and  (52),  which  then  become 


X 


(53) 
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and 


X 


(54) 


To  obtain  the  probability  distribution  for  x,  we  assume  that  all  initial 
phase  angles  dj,  -{n  +  0^)  ^  rr  •_  {n  -  are  equally  probable;  that  is,  the 

1 

probability  distribution  for  0-  is  uniform  and  of  value  p  (dj)  =  -  (see 

In 

Fig.  9). 

Since  the  use  of  o  as  the  dependent  va l  iable  instead  of  d-  is  simply 
a  linear  shift  of  the  angular  coordinate,  the  probability  distribution  for  8  itself 

1 

is  uniform  and  of  value  p  (<S)  -  —  . 

2n 

The  probability  distribution  p  (\)  consists  of  two  parts,  one  discrete  and 
the  other  continuous.  The  discrete  part  consists  of  an  impulse  at  X  =  0 
magnitude  ~  .  Tliis  can  be  seen  hy  considering  Fig.  9,  which  shows 

that  initial  phase  angles  d^  sucli  that  n  ^  ...  ,  ^  f  already  lie  within 

the  prescribed  distance  of  and  are  therelore  regarded  as  synchronous 
and  requiring  zci  o  time  to  pull-in. 

The  continuous  part  of  the  proliability  distribution  is  determined  below. 

It  is  again  convenient  to  break  this  problem  into  two  sections  according  to 
the  two  possible  ranges  of  dj  . 


-(77  .  d^^)  ^  dj  {0^^  -  (  ) 


(e,  .  0  _ 


(.7  -  d 


O 


) 


Fig.  9,  Probability  Distribution  for  Initial  Phase  Angle  0. 


PROBABILITY  DISTRIBUTION  p  (X)  FOR  ,  -(n  +  6^)  <  (9^  <  {6^  -  () 
Using  Eq.  (53) 


or 


B  1  /  S  \ 

X  =  —  +  —  InIA  cot  —  +  a  I 
A  A  \  2  ) 


(55) 


(56) 


38 


where 


B  = 


is  a  constant  dependent  on  the  mistuning  and  e ,  and  does  not  contain  the 
variable  S. 

Since  X  is  a  continuous  single-valued  function  of  0-  for  the  range  of 
values  of  0j  considered,  the  probability  that  the  variable  X  lies  between 
the  values  x  and  x  +  dx  is  equal  to  the  probability  that  the  variable  0j  lies 
between  the  values  0  and  0  +  dO  .  This  is  also  equal  to  the  probability  that 
the  variable  S  lies  between  the  values  0  and  6  +  dO  as  previously  explained; 


p  (x  <  X  <  X  +  dx)  =  p  (0  <  0-^  <  0  +  dO) 


=  p(O<S<0  +  dO) 


In  differential  notation 


p  (X)  dX  =  +  p  (5)  dS 


or 


1 

p  (X)  =  +  p  (5)  . 

dx 

d8 


(57) 
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Since 


P  (3)  =  P  (0i)  = 

‘  2  7T 


we  require  only  —  in  order  to  evaluate  p  (X)  using  Eq.  (57). 
dS 

Using  Eq.  (56) 


We  now  require  the  trigonometric  function  in  Eq.  (59)  in  terms  of  X  . 
Using  Eq.  (58) 


AX  -  B  = 


40 


or 


Now 


.AX-B 


«AX-B 

e  —  a 


/  5  \ 

=  I  A  cot  —  4-  a  I 

\  2  } 


=  cot 


cosec  —  =  1-1-  cot  — 
2  2 


(60) 


cosec*^  —  =  1  + 

2 


C-^) 


(61) 


Using  Eqs.  (60)  and  (61)  and  substituting  in  Eq.  (59) 


dX 

Us 


1  + 


.AX-B 


(62) 


dX 

Substituting  in  Eq.  (57)  for  —  obtained  in  Eq.  (62)  and  recalling 

d5 


that  p  (S)  =  — 

277 


p  (X)  =  - 
n 


:.AX-B 


(^) 


1  + 


(63) 
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PROBABILITY  DISTRIBUTION  p  (X)  FOR  ((9^  +  0  <  -  0^) 

Using  Eq.  (54) 


X 


(64) 


or 


1  /  S  \ 

X  =  —  +  —  In  (A  cot - a  I 

A  A  \  2  / 


(65) 


where 


-  In 


and  is  a  constant  dependent  on  the  mistiming  and  e ,  and  does  not  contain  the 
variable  S.  The  similarity  between  Eqs.  (65)  and  (66)  is  apparent,  so  to  obtain 
p(X)  for  +  ( <9-^<(77  -  9^)  we  just  change  B  to  B^  and  -a  to  +a  in 
Eq.  (63)  and  we  have  the  probability  distribution 

p(X)  =  — 
n 


.AX-B 


1 


1  + 


(66) 
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FINAL  EXPRESSION  FOR  THE  PROBABILITY  DISTRIBUTION 


Using  the  results  obtained  from  Eqs.  (63)  and  (66)  and  recalling  that 
p(X  =  0)  =  -1 

77 


P 


(X)  =  -  for  X 


0 


(67) 


and 


p  (X) 


gAX-B^ 


for  X  >  0  where 


B  =  In 


(69) 


and 


B 


1 


In 


cot 


( 

1 


(70) 
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The  expression  for  the  probability  density  function  given  by  Eq.  (68) 
may  now  be  rewritten  with  0^  as  the  dependent  variable  by  recalling  that 


a  =  sin0^  and  A  =  cos0^ 


Equation  (68)  reduces  to 


p  (X )  = 


2  a 

cos 


1 


2  TT 


cosh  (AX-B)  - 


1 

-f-  - 

cosh  (AX— B^)  +  sin  6^ 


(71) 


where 


sin  — 
2 


and 


B 


1 


(73) 
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Since  later  appendixes  are  concerned  with  the  mean  and  the  accumulative 
distribution,  it  is  necessary  to  verify  that  the  integral  of  the  probability  density 
does  in  fact  equal  unity,  that  is 


p(X)dX  =  1 


Using  Eqs.  (67)  and  (68),  if 


oo 


i 


p(X)  =  1,  it  is  required  to  show  that 


I 


oo 


1  + 


^AX-B 


+ 


AX-B 

e 


dX 


TT  —  € 


(74) 


To  evaluate  the  first  term  in  Eq.  (74);  i.  e. , 

eAX-B 

-  dX 

eAX-B  _  ,y 

we  use  the  substitution 


u 
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Then 


Similarly  the  second  term  in  Eq.  (74);  i.  e. , 


f 


.AX-B^  dX 


—  tan 


l(e ®  +a 


1  + 


/  AX-B^ 
re  +  a 


(76) 


Using  Eqs.  (75)  and  (76),  the  left-hand  side  of  Eq.  (74)  is  equal  to 
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To  show  that 


P  (X)  =  1,  it  remains  to  prove  that 


1 


oo 


can 


-1 


(77) 


Recalling  that 


e 


-B 


€ 

cot  —  +  a 
2 


) 


-1 


(from  Eq.  (56)) 


and 


e 


(from  Eq.  (65)) 


Equation  (77)  becomes 
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For  conveniejice  Let 


s  =  cot  — 

2 


(79) 


u  =  tan 


V  =  tan 


.1 

y  As  +  a  J 
-1  /a  +  as  \ 

U®  -  V 


(80) 


(81) 


Taking  the  tangent  of  each  side  of  Eq.  (78)  we  get 


tan(u  +  v)  =  tan  e 


(82) 


Now 


tan(u  +  v)  = 


tan  u  -I-  tan  v 
1  —  tan  u  tan  v 


(83) 


Using  Eqs.  (80)  and  (81) 


A  —  as  A  -1-  as 
-1- 


tan  (u  +  v)  = 


As  +  a 


As  —  a 


/A  —  as\  /A  +  SLi 
I  As  +  a  y  y  As  — 


(A  —  as)  (As  —  a)  +  (A  +  as)  (As  +  a) 
(As  +  a)  (As  -r  a)  —  (A  —  as)  (A  +  as) 


(84) 


2s  (A^  +  a^) 

(s^  -  1)  (a2  +  a^) 


2s 


(s^  -  1) 
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Since  s 


=  cot  —  (from  Eq.  (79),  we  have 


tan  (u  +  v) 


2  cot  — 
2 


tan  (u  +  v) 


2 


1  - 


tan  € 


which  completes  the  proof  required  by  Eq.  (82). 


APPENDIX  m 


DETERMINATION  OF  MEAN  FOR  PULL-IN 
TIME  T  OF  A  FIRST-ORDER  PHASE-LOCKED  LOOP 


Appendix  II  shows  that  the  probability  distribution  for  X  is  given  by 


P(X) 


—  for  X  =  0 
n 


(85) 


and 


P(X) 


1 

77 


gAX-B 


+ 


AX-B^ 

e 


(86) 


for  X  >  0 , 
where 

X  =  KT 

K  =  open-loop  gain  of  the  loop  (a  constant) 

T  =  pull-in  time  of  the  loop 

The  impulse  at  X  =  0  does  not  contribute  in  any  way  to  the  mean;  hence 
it  can  be  ignored.  The  mean  of  X  is  therefore  given  by 


X 


X  p  (X)  dx 


(87) 
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Using  Eq.  (86)  gives 


X  = 


1 

77 


X 


X 


dX 


(88) 


It  is  not  obvious  from  Eq.  (88)  which  of  the  many  possible  substitutions 
will  give  a  convenient  solution.  The  trick  in  obtaining  the  solution  to  Eq.  (88) 
is  to  recall  how  the  expression  for  p  (X)  was  derived  in  the  first  place. 

Equations  (55)  and  (64)  were  used,  originally,  to  determine  p  (X);  i.  e. , 


If  we  substitute  X,  as  given  in  Eq.  (89),  in  the  first  term  of  Eq.  (88), 
and  X,  as  given  in  Eq.  (90),  for  the  second  term  in  Eq.  (88),  then  Eq.  (88) 
reduces  to 
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where 


(91) 


Consider  the  first  integral  in  Eq.  (91);  i.e. , 


In 


6 

cot  - 


2 


(tt  -  26^  -  6)  + 


In 


cot 


d8 
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Consider  the  second  integral  in  Eq.  (91) 


n  +  20  —  f 
o 


dS 


77  +  20  —  €  ^ 

O 


77  +  Zi 

/ 

•'5=0 


In 


/  ^  \  ^  \ 

I A  cot  —  +  a  1  +  In  I A  cot  _  +  a  I 


d5 


(93) 


=  In 


(A  cot  —  +  a  J 

2  I 


7T+26^-e 


(77  +  20^  -  0  + 


In  (  A  cot  —  +  a  ] 
5  =  0  ^  ^ 


d5 


Using  Eqs.  (92)  and  (93)  and  substituting  in  Eq.  (91) 


X  = 


1 

2  77  A 


/ 

I  A  cot  — 
\  2 


(tt  -  20^  -  f)  +  In 


(A  cot  —  +  a  1 

2  j 


(77  +  20^ 


-0 


+ 


c 


(  ^ 

In  I A  cot  — 

V  2 


77  +  2l 

a j  dS  +  I 

’  h  =  o 


7T+26^-e 


(94) 


In 


IA  cot 


5  \ 

—  +  a  1  d5 


Since  A  =  cos  6^ ,  and  a  =  sin  6^ 


S 

A  cot  —  —  a  = 
2 


sin 


5 

T 


(95) 
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and 


8 

A  cot - h  a  = 

2 


(96) 


Using  these  results  gives  for  the  second  part  of  Eq.  (94) 


which,  after  considerable  rearrangement  of  limits,  reduces  to 


1 

2  77  A 


2 


€ 


In  cos  (A)  dA  —  4 


In  sin  A  dA 


(97) 
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Equation  (94)  now  becomes 


If  f 

expression 


is  small,  the  last  two  terms  in 
for  X  reduces  to 


Eq. 


(98)  may  be  neglected  and  the 


(99) 
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APPENDIX  IV 


CUMULATIVE  DISTRIBUTION  FOR  THE  PULL-IN  TIME  T 

The  probability  density  for  the  pull-in  time  T  is  given  by  Eqs.  (67)  and 
(68);  i.  e.  , 


p  (X)  =  1  for  X  =  0 

n 


(100) 


and 


P 


gAX-B 


+ 


gAX-B^ 


(101) 


for  X  >  0 , 
where 


X  =  KT 


and 


K  =  open- loop  gain  of  the  loop. 
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The  cumulative  distribution  for  X  is  given  by  the  probability  that  X  is 


less  than  some  specified  value,  say  X^  ;  i.  e. , 


P  (X  <  X^) 


p  (X)  dX 


Using  Eqs.  (100)  and  (101) 


P  (X  <  X^) 


Consider  the  first  integral  in  Eq.  (102) 


I 


1  - 


gAX-B 


Let 


^AX-B 


a 


A 


u 


h 


du 


1  +  u 


2 


(102) 
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where 


=  (^) 


and 


ax_-b 


U2  = 


=  I  tan  ^  u  I  =  tan 


“2 


,-l| 


e-B_a 


(103) 


Using  a  similar  substitution  the  second  integral  in  Eq.  (102)  reduces  to 


I 


2  ■ 


(104) 


Using  the  results  of  Eqs,  (103)  and  (104)  we  obtain 


58 


In  Appendix  II  (Eq.  (77)  onward),  it  is  shown  that  the  last  two  terms  in 
Eq.  (105)  were  equal  to  f ,  so  that  Eq.  (105)  reduces  to 


p  (X  <  X^)  = 


tan 


-1 


+  tan 


-1 


AX  -B^ 
e  °  +  a 


(106) 
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APPENDIX  V 


A  NOTE  ON  INSTANTANEOUS  PHASE  AND  INSTANTANEOUS  FREQUENCY 

This  Appendix  explains  the  terms  "instantaneous  phase  difference"  and 
"instantaneous  frequency  difference. " 

Suppose  we  are  given  a  sinusoidal  waveform 

El  =  sin(9^  (t)  (107) 

We  define  its  instantaneous  phase  as  (t)  and  its  instantaneous  frequency 
as  the  time  derivative  of  its  phase  angle 


(t) 


(t) 

~dt 


(108) 


For  example,  consider  the  case  where 


(0  =  t  +  a  (109) 

which  corresponds  to  a  sinusoid  with  constant  angular  frequency  and  an  initial 
starting  phase  at  time  t  =  0  of  a. 

Equation  (109)  gives  its  instantaneous  phase  angle  as  a  function  of  time 
and  the  time  derivative  of  Eq.  (109)  gives  its  instantaneous  frequency  as  a 
function  of  time. 


di  (t)  = 


(110) 
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For  this  case  the  word  instantaneous  is  redundant,  since  we  are  considering 
a  sine  wave  of  constant  frequency. 

Suppose  now  we  are  given  two  sinusoidal  waveforms 


=  sinSj  (t) 

(111) 

=  sin02  (t) 

(112) 

We  define  their  instantaneous  phase  difference  as 


d  it)  =  01  (t)  -  02  (0 


(113) 


and  their  instantaneous  frequency  difference  as 


0  (t)  =  01  (t)  -  02  (t) 


(114) 


For  example,  consider  the  case  where 


01  (t)  =  (Ui  t  +  a 


and 


02  (t)  ~  (o^  t  +  ^ 


corresponding  to  two  sinusoids  with  different  angular  frequencies  and  different 
initial  phase  angles. 
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Using  Eq.  (113),  their  instantaneous  phase  difference  is 


9  (t)  =  (t)  -  02  (0 

=  («Uj  t  +  a)  —  (ct)2t  +  /8)  (^^5) 

=  (cjj  -  (U2)t  +  (a  -  /8) 

Using  Eq.  (114),  their  instantaneous  frequency  difference  is 


h  (t)  =  0j  (t)  -  02  (t) 


(116) 


=  (wj  -  W2) 


The  definitions  of  instantaneous  phase  difference  and  instantaneous 
frequency  difference  given  by  Eqs.  (113)  and  (114)  are  always  valid  and 
independent  of  the  actual  form  of  0j  (t)  or  02  (t)  . 
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